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Abstract 



I We consider the situation of a two-level quantum system undergoing a continuous 

Qh' indirect measurement, giving rise to so-called "quantum trajectories". We first describe 

fin I these quantum trajectories in a physically realistic discrete-time setup and we then justify, 

by going to the continuous-time limit, the "stochastic Schrodinger equation" attached to 
this model. We prove return to equilibrium properties for these equations. 



1 Introduction 

> , 

, Recent experiments of continuous measurement in quantum mechanics (Haroche's team in 

particular), or more precisely in quantum optics, have put into evidence the random evolution 
^f) • of the state of a quantum open system [HarH IHar2] . In particular, one has experimentally 

. observed "quantum jumps". These experiments allow to study the evolution of a quantum 

I system interacting with some environment. They are based of the principle of indirect mea- 

surement on the environment, in order not to perturb the evolution of the small system 
[BaHlBi3l[Dm[Di2llWMTllWM2llB?Pj . 

The stochastic models attached to these phenomenons are described by stochastic dif- 
^ I ferential equations, called "Stochastic Schrodinger Equations" or also "Belavkin Equations" 

^ '• [BiGl [BaBl iBiTl [Ba2l [Bi3l iBiZl [BiHl iBeMl [BeMl [ED EUl Eia IWMTl IWM2llB?P] . Their 

solutions are called "quantum trajectories", they describe the evolution of the state of the 
small open quantum system. The stochastic differential equations which are usually obtained 
in this context are of two different types. 

Either they are of "jump-type" : 

dpt = L{pt)dt + (tM^^ - pt] [dNt - TT[J{pt)\dt] . (1) 



Tr[J{pt)] 

where {Nt) is a stochastic counting process with stochatic intensity /g Tr[j7'(/9s)]ds. The op- 
erator L corresponds to a Lindblad type operator and the operator J describes the evolution 
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of the system during the quantum jumps. This equation describes experiments which are 
called "direct photon detection" (observation of the photon emission by an atom excited by 
a laser). 

Or it can be an equation of diffusive type: 

dpt = L{pt) dt + (Cpt + ptC* - Tv[{C + C*)pt]pt) dWt, (2) 

where (Wt) is a standard Brownian motion. In quantum optics, this equation describes 
experiments called "Heterodyne or Homodyne detection". 

In the usual literature, obtaining and justifying rigorously these equations makes use of 
Quantum Filtering Theory |Ba3[ IBaH IBeMl IBGMj . It is the quantum probability version of 
the usual filtering technics, it makes use of fine Quantum Stochastic Calculus and heavy von 
Neumann Algebra Theory. Others approaches are based on classical probability and use of 
Instrumental Process and notion of "a Posteriori State", cf |BaGl [Ba2t IMoRl IBaH] . 

A maybe more intuitive and more physical approach for these equations is to start from 
a discrete-time procedure, that is, repeated quantum interactions with measurement of the 
environment ( [APlj ). One can then obtain the stochastic Schrodinger equations by passing 
to the limit to a continuous-time model |PeH IPe2| . 

In this article, we come back and apply results obtained in [Pelj and [Pe2) . in which 
Belavkin equations are obtained with this approach. Here, we obtain the description of the 
stochastic Schrodinger equations for a two level system in contact with a spin chain. We 
adapt the result of [Pel] and [Pe2j in order to describe the quantum trajectories in terms 
of wave function (in [Pel] and |Pe2] . the stochastic equations for the evolution of density 
matrices have been derived from the approximation procedure). Next, for a special model, 
we show a property of return to equilibrium of the solution. 

2 Discrete-Time Quantum Trajectories 

In this section we describe the physical model and the mathematical setup of indirect repeated 
quantum measurements. We describe the evolution of the small system undergoing successive 
measurements through the "discrete quantum trajectories". 

2.1 Repeated Quantum Interactions 

The physical situation is the following. A quantum system, with state space Us (often called 
small system for it is in general finite-dimensional and small compared to the environment) is 
undergoing repeated interactions with a chain of quantum systems (S'n*'^- This is to say that 
we consider an environment which is made up of a sequence of identical copies of a quantum 
system, each with state space T-L. Each piece Ti of the environment is going to interact, one 
after the other, with the small system Us- This interaction lasts for a time duration r and is 
driven by a total Hamiltonian Htot on T-Ls T-L- Hence, each interaction is described by the 
unitary operator 

U = e~'^"'°' 

on Tis ®T-i- In the Schrodinger picture, if p denotes any initial state on the tensor product 
T-Ls ® Ti then the evolution of the state after this interaction is given by: 

p^UpU* . 
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After this interaction, the systems T-Ls and % stop interacting together, the system T-Lg comes 
to meet a second copy of 7i and they interact together in the same way as before (that is, with 
following the same unitary operator U). And so on ... the small system interacts repeatedly 
with each of the independent copies of %. 

Let us develop the mathematical framework which allows to describe these repeated quan- 
tum interactions. We follow the setup of the article [APlj . in which these models and their 
continuous limit were first introduced. 

The state space describing the whole game is 

r = ^5 ® (g) ^fc , (3) 

where each is a copy of the Hilbert space %. We have to be clear about what the above 
countable tensor product means: 

r$ = (g) T^fc . 

km* 

Recall that a countable tensor product of Hilbert spaces can only be defined with respect to 
a choice of a particular unit vector Uk in each copy %k (the so-called the stabilizing sequence 
of the countable tensor product). In our case, we assume that T-Lk is finite dimensional and 
we choose an orthonormal basis 

{X'- i G A^U{0}} 

where is a set of the form {!,... A^}, which is the same for each T-Lk. A particular role is 
played by the vector which has to be considered as a reference vector for the system 
as we choose the stabilizing sequence to be Uk = X^ for all k. 

Denote by X^, the basis vector X* but leaving in the A;-th copy T-Lk ofH. Then an Hilbertian 
orthonormal basis of T$ is given by all the tensor products ^k^k where all the vectors Vk are 
equal to X^, except for a finite number of them which might be equal to some AT^*, ik G AA. 
This stands for a definition of the countable tensor product = <Sik&n*T^k- 

The repeated quantum interaction setup is based on two elements: the time length r and 
the Hamiltonian Htot which describes each basic interaction. Consider the unitary operator 
U = exp(— iri^tot) acting on T-Ls H and consider the unitary operator Uk on T which acts 
as U on Tis (8) Tik and which acts like the identity operator on the other copies Tik'- This 
operator Uk describes the effect of the A;-th interaction. 

The unitary operator 

Vk = Uk...Ui 

describes the effect of the k first interactions. Indeed, if p is any initial state on T, then 

VkpVk* 

is the state of the whole system (small system -|- environment) after k interactions. 
Define the elementary operators a*, i,j £ J\f Ci {0} on T-l by 

aix'' = 5i^kX^ . 

It is useful for further computations to notice that in Dirac notation a*- = \X^){X^\. We 
denote by aUn) their natural ampliation to T<I> acting on the n-th copy of Ti only. 
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Clearly, U can always be written as 

u= 

ijeA/'u{o} 

for some operators Uj on "Hq such that: 

kej\fu{0} feeA/'u{o} 

With this representation for U, it is clear that the operator representing the n-th inter- 
action, is given by 

E U'j^aiin). 

i,jeAfu{o} 

With these notations, the sequence (Vn) of unitary operators describing the n first repeated 
interactions can be represented as follows: 

Vn+l = Un+1 Vn 

= E U'j0a]{n + l)Vn. 
i,jeA/'u{o} 

But, inductively, the operator Vn acts only on the n first sites of the chain T^, whereas the 
operators a* (n + 1) act on the (n + l)-th site only. Hence they commute. In the following, 
we shall drop the (X" symbols, identifying operators like a* (n + 1) with ly^^ a*(n + 1), the 
operator f/j with Uj It^, etc. This gives finally 

Vn+i= E U'jVnaiin + 1). (4) 

ijGA^U{0} 

On T^, one vector plays a particular role, the vector 

n = 0kXt . 

For any bounded operator K on F, we define the operator Eo[-ftr] on Tis as the unique operator 
on Tis such that, for all trace-class operator p on Tis we have 

Tr^g(pEo[K]) =Trr((p® |0 >< n\) K) . 

That is, Eo[-f^] is the partial trace of K with respect to the state |0 >< on T^. 

We then have the following fundamental action of the repeated interactions, when re- 
stricted to the small system. 

Theorem 1 (cf |APlj ) The effect of the repeated interaction dynamics when restricted to 
Us is given as follows. For all observable X on lis, for all n G N, we have 

MV:{X®I)Vn\ = L^{X), 

where L is a completely positive map on Tis whose Krauss decomposition is 

L{X) = Y,{Uf)*XU^. 

Any (discrete) semigroup (L") of completely positive maps can be obtained this way. 
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Note that the completely positive map L defined above acts on observables. It also induces 
a completely positive "dual map" L* acting on states as follows: 

L*{p) = Y,Ufp{Uf)* (5) 

and which satisfies 

Ti{pL{X)) = T,{L*{p)X) 

for all state p and all bounded operator X on 'Hg- Recall the usual notion of partial trace 
defined as follows. 



Definition-Theorem 1 Given any state a on a tensor product T-L® K, then there exists a 
unique state rj onlH. which is characterized by the property: 

Tr[r?X] =Tr[a {X I)] , 

for all X G B{'H). The state rj is denoted by Trx;(a) and is called the partial trace of r/ with 
respect to IC. 

With these notations we have the following result. 

Theorem 2 For every state p on T-Ls md all n E M we have 

ivt$(f„(p® |o X o|)v;*) = (L*f(p) . 

Proof: We have, for all X bounded operator on Hs, 

Tr((L*r(p)X)=Tr(pL"(X)) 

= IV(pEo[K(X®/)V;]) 

= Tr ((p ® \n >< n\) v*{x ® i)Vn) 
= Tr {Vnip ® \n >< n\)v* {X (g) J)) 
= IV (ivt* {Vn{p ® \n >< n\)v*) x) . 

This proves the announced result. □ 



2.2 Repeated Quantum Measurements 

We now somehow consider a more complicated procedure. After each interaction is finished, 
the piece Hk of environment which has just finished to interact with T-Ls is undergoing a quan- 
tum measurement of one of its observables. The random result of this quantum measurement 
will give some information on the state of the whole system and in particular on the state of 
Hs- The so-called quantum trajectory is the random process we obtain this way, by looking 
at the knowledge we have of the state of Hs after each measurement. 
Let A be any observable on H, with spectral decomposition 

p 
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the Aj's being the eigenvalues, the Pj's being the eigenprojectors. We consider the natural 
ampliations of A which defines an observable on T by making A acting on the A;-th site T-L^ 
only: 

fc-1 

A^ = ^I(^A(^ (g) / 

j=0 j>k+l 




with obvious notations. 

As a consequence, if p is the state of T then a quantum measurement of the observable 
A^ gives the values Xj with probability: 

P[to observe A^] = rr[pPj= ], j e{l,... ,p} . 

If we have observed the eigenvalue Xj for the observable A'', the new state of the system is 

pfc^pfc 
~ Tr[pPf] ' 

This principle is the so-called "von Neumann projection postulate" . Now, if we perform an- 
other measurement of the observable A^ we obtain P[to observe Xj] = 1. As a consequence, 
a naive repeated measurement operation gives no information on the evolution of the sys- 
tem. The repeated measurement procedure has to be combined with the repeated interaction 
procedure in order to give non-trivial informations on the behavior of the system. 

The quantum repeated measurement principle is the combination of the measurement 
principle and the repeated quantum interactions. Physically, this means that each copy T-Lk 
of T-L interacts with T-Ls and we perform a measurement of A^ on T-L^ after it has interacted 
with T-Ls- After each measurement we have a new (random) state of the whole system, given 
by the projection postulate. This is the so-called discrete quantum trajectory. 

More precisely, the initial state on T is chosen to be of the form 

where p is any state on T-Lq and each rji = rj is a reference state on Ti. We denote by fi^ the 
state representing the new state after the k first interactions, that is. 

Let us now define the probabilistic framework in order to describe the effect of the successive 
measurements. We put $7 = {1, . . . ,p} and on Cl^ we define the cylinders of size k: 
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Wc endow 17 with tlie u-algebra T generated by all these sets, this is the cylinder a-algebra. 
Note that for all j, the unitary operator Uj commutes with all the projectors Pj' such that 
k ^ j. Hence, the state of the system after k interactions and k measurements which have 
given the respective values Ajj , . . . , Aj^. is (up to normalization by the trace) 



= pt^...plUk... u,^i{u,r ... {UkTPl ---Pi 

~ Pik ■ ■ ■ ■ ■ ■ Pik ■ 

We denote by fi{ii, . . . ,ik) the quantity 



Pt...Plf^kPl,...Pt- 

By the Kolmogorov Consistency Theorem we can define a probability measure P on (fi^ , J-) 
only by specifying 

P[Ai,,...,iJ =TV[/i(ii,...,ifc)]. 
We also define a random sequence of states on T by 

^ I — > Pk{!^\...ojk)- 



Tr[/i(a;i . . .Wfc)] 



This random sequence of states is our discrete quantum trajectory and the operator . . . , i^) 

represents the state of the system, after having observed the results Ajj , . . . , Aj^ for the k first 
measurements. This fact is made precise in the following proposition. 

Proposition 1 Let (pk) be the above random sequence of states we have, for all uj G ^ 

p^+^\Uk+iPkHiUk+irp^+^\ 



Tr 



This proposition is obvious but summarizes the quantum repeated measurement principle. 
The sequence pk is the quantTim trajectory, showing up the effect of the successive measure- 
ments on r. The following theorem is an easy consequence of the previous proposition. 

Theorem 3 The sequence {p^)n is a Markov chain, valued in the set of states ofV. It is 
described as follows: 

// p" = 6n then p""*"^ takes one of the values: 

Pi^^Un+l On {Un+l)* Pi^^ 



T,[Un+ien{Un+irP^+^] ' 

with probability Tr [C/„+i On (C/n+i)* Pi""^ ] ■ 
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The most interesting behavior of the Markov chain of states above is obtained when one 
restricts it to the smah system Tig- This way we obtain a quantum trajectory on the states 
of Us by considering the sequence of random states on Us'. 

Pn{uj) = TrT$(pn(w)) . (6) 

This defines a sequence of state on T-Ls which contains the "partial" information given by the 
measurement and we have the fohowing theorem which completely describes the behavior of 
this random sequence. 

Theorem 4 The random sequence defined by formula ^ is a Markov chain with values in 
the set of states on T-Ls- If Pn = Xn then Pn+i takes one of the values: 

Tin [{I ® Pi) UjXn ^ r,)U* {I ® Pi)] .^^ 
Ti[U{xn®r^)U*{I®P,)] ' ' 

with probability TV [U{xn r])U* (I <^ Pi)] ■ 
The expectation of pn satisfies 

E[p„] = (r)"(/,o), 

where L* is the completely positive map described in Theorem\^ 

Proof: Assume, by induction, that pn is given. This means that TiT^iPn) = Pn- The 
next step of the quantum measurement gives (by Theorem [3]) 



Pn+l 



P^+^Un+lPn iUn+l)*Pj'^' 



Tv[Un+l9„iUn+l)*PP'-'] ' 

for some i. Hence, we have to compute 

TlMPP+'Un+l Pn (C/„+l)*/^"+l) . 

Decomposing, with obvious notations, the space T into Hs T<I>jo,n] T~in+i ^^[n+2,+oo[) 
one notes that, by induction, the state pn is of the form 



fc>n+2 



where On is a state on Us ® T^[o,n]i satisfying 



TrT$[o.„] (^n) = Pn 



Hence, for all X, bounded operator on Us-, we have 

Tl{TTMPj^^'Un+lPn{Un+irP:^')X) = 

= Tr {{P^+^Un+l Pn {Un+l)*PP''') {X ® I)) 

= Tr {Un+l Pn {Un+lY {X ® /[o,„] ® Pj'^^ ® /[„+2,+oo[)) 



Tv\Un+l \en®r]® r]\ {Un+lf {X ® I^^^n]® P'i 

k>n+2 



.71+1 



71+2, +oo[ 

■n+l\ 



TV (( 0„ » r/) {Un+lT{X I[o,„] PP+^)Un+l) . (7) 
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But Un+i acts only on Hs (8) Tin+i, hence the operator {Un+i)*{X I[o,n] Pl^~^^)Un+i is 
an operator on Tis T-ln+i ^^[o,n] (note the interchange of space, for simphcity of the 
notations) which is of the form 

{{Un+inXCSPr^')Un+l)CSllo,n]- 

Hence, the quantity ([7]) is equal to 

TV (TrT<i.[o,„, (e„ ® r?) i^"+i)C/„+i) . 

But TrpjO/re] {Gn®ri) is equal to Trr|0,n] ( ^n) = Pn®V ■ This gives finally 

Tr ( TrT*(i^"+^C/„+i p„ (C/n+i)*/^^!) = 

= TV {{Pj'+^Un+i {pn v) X) . 

But in this expression, the index n + 1 plays no more role and the expression above may as 
well be written 

TT{{PiU{pn^V){UTP,)X) 

on Tis T-L. This proves the first part of the theorem. 

Let us check, the one concerning the expectation of pn- Note that the expectation of pi 
is equal to 

p 

= TrniUipo CS> r])U*P,Pi) for P, acts on U only 
1=1 

p 

= Trn{U{po(^v)U*YPi) 
1=1 

= TTn{U{po0v)U*) 
= L*{po). 

We conclude easily by induction. □ 

Thanks to the above description we can express a discrete-time evolution equation for the 
quantum trajectories. Let us put 

C^{p) = Eo [(I ® Pi) U{p ® 7])U* {I Pi)] , 

i = 1, . . . ,p. We then have for all uj G and all A; > 0: 



i=o Tr[£i(pfe)(tj)] 



where if (w) = hiiVk)- 
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2.3 The two-level atom model 



In this section we specialise to the case where Hq = C^, this is the so-called two-level atom 
model. In most of the physical applications that we have in mind, the interacting system is 
also of the form % = We denote by Xq , Xi an orthonormal basis where the reference 
state r] is diagonal: 

' Vo 
VI 

Let X be any orthonormal basis of Hq. For describing the interactions between T-Lq and 
Ti we choose (gi Xq,X Xq, fl (8> Xi,X (8> Xi as an orthonormal basis of T-Lq In such 

a basis, the unitary operator U , describing the elementary interaction, can be written as a 
2x2 matrix with coefficients being operators on T-Lq. That is, we can write U as: 



ri 



U 



Let A be an observable of % on which we want to perform a measurement. It can be written 
as ^ = Aq-Pq + AiPi where Aj are its eigenvalues and Pi the corresponding eigenprojectors. 
Let {PI i)k,i=o,i be the matrix elements of the projector in the basis Xo,Xi. Put 



k,l=0,l 



Then, if denotes the state of the system Tis after the k-th measurement, the state Pk+i 
takes one of the two possibles values 



We denote Pk+i = Tr[£o(/Ofc)] or qk+i = TT:[Ci{pk)] the corresponding transitions probabilities. 

In the rest of the paper, we concentrate on a special case of environment, where rj = 
\Xq){Xq\. This situation corresponds to a model of heat bath at zero temperature, see [AtPej 
for more explanations and for positive temperature models (let us just stress that this choice 
is crucial and that positive temperature gives rise to completely different continuous-time 
behaviours). In this situation, the discrete quantum trajectory can be described in terms 
of pure states. More precisely, if the initial state of Hq is pure, the random sequence (pk) 
remains pure. This way, we can describe the evolution of Hq with a random sequence of 
vectors (wave functions). 

Proposition 2 Let Tio = H = 'C? and rj = \Xq){Xq\. Let (pi^) be the discrete quantum 
trajectories corresponding to the indirect measurement of an obervable A. 

If Po is a pure state, that is po = \ipo){'ipo\, oind if the measurement is non-trivial (A is 
not a multiple of the identity), then the state of the small system pn is always a pure state. 
In other terms, there exists a random sequence of wave functions (|V'n)) such that \\tpn\\ = 1 
and such that pn = \ipn) {i'nl , for all n G N. 

The sequence {\ipn)) is also called a discrete quantum trajectory. 

Proof: Since we work in 2-dimension and since A is not a multiple of identity, we have 
A = XqPq -\- XiPi where the Pj's are one dimensional projectors. Thus there exist two vectors 
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Oi, i = 0,1 such that Pi = \ai){ai\. Now, let po = |V'o)('0o|) after the first measurement if 
we have observed the eignevalue Aj, the non normahzed state describing the experiment is 
described by 



Eo[l0\ai){ai\ U{\^Po){H<^\Xo){Xo\)U* I0\ai){ai\\ 

,0 \l^\ai){ai\ (ui{\i;o){'ilJo\{Ky ^ ai\Xo){Xo\al^ I®\ai){ai 



k,l u,v 



k.v 



k,v 

E 



^{ai,X^)U',i;o) (Y.{a^,X^)Uyo 

k / \ V 



fX" \ai){ai 



Now, by normahzing the vector Ylki^i^ ■X^)^k'^o^ it is straightforward that we get a vector 
V'l such that pi = I- Next, by induction we can construct a sequence ■0„ such that 

Pn = \'4^n){'^n\ for aU n. □ 

Remark: Such a property is at the basis of the use of "quantum trajectory theory" 
for numerical simulations of Lindblad master equations. Numerically, the description in 
terms of pure states reduces the number of parameters to control (in comparaison with den- 
sity matrices) . We recover the "deterministic" dynamic by taking the expectation, that is, 
^[|V'ril(V'n|] = ^"'{po)- In the continuous time version, similar properties are called "unrav- 
ellingj' of master equations and simulations make use of technics called "Quantum Monte 
Carlo simulations". 

Now we wish to obtain discrete-time evolution equation which describes the stochastic 
evolutions of discrete quantum trajectories. To this end, let us introduce some notations. Let 
Pq be the projector on uq = (/i, z^), with ||ao|| = 1 and Pi the projector on ai = {u, —fl). Let 
define the following functions acting on vectors 



Mm 
Mm 



\[^^U^ + lyU'^)] ^P) 



Then, the dynamic of (V'n) can be described by the equation 



\mM^)))\\ 



,k+i, 



CO) + 



l-^i(I^.M))ll 



,k+l 



(9) 



for all uj G S^. This equation corresponds to equation ([8]) for a two level system in terms of 
wave functions, i.e. the sequence (|V'fc)(V'fe|) satisfies equation ([8]). 

In the next section, we will describe the continuous time version of these equations. To 
this end, we aim at considering this discrete-time model but depending on a time-length 



Unravelling means the description of a wave function stochastic process (V't) such that E[ji/'t)(V't|] = e'-^(po) 
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parameter r which we shah make tend to 0. That is, we want to pass from a discrete time 
interaction model to a continuous time one. This way, we shall obtain the classical Belavkin 
equations for quantum trajectories associated to continuous measurement. In the litterature, 
these equations describe a model where a two-level atom is in contact with a photon-stream. 

Let = :^ be the time of interaction between the small system and one element of the 
environment. Let us denote by U{n) the unitary operator associated to each interaction, it 
now depends of the time of interaction. If we had no measurement process on the environment, 
we will be back to the problem of going from a discrete-time repeated quantum interaction 
model, to a continuous time one. This problem has been completely studied in |AP1| . In 
their article they show that, in order to get a limit evolution when r goes to 0, we have to 
ask the operator U{n) to satisfy certain renormalization conditions. They have shown that 
the coefficients Uj{n) must follow well-defined time scaling in order to obtain a non-trivial 
limit. Namely they have shown that the operators Vint] — U{[nt]) . . . Ui, t > 0, converges to 
an evolution {Vt)t which is a continuous operator process. This process naturally satisfies a 
quantum Langevin equation which represents the evolution equation of the small system -|- 
bath. 

Our continuous measurement procedure does not differ much from their approach, except 
that we perform a measurement on the environment after each interaction. This is why we 
have to keep the same normalization for the coefficients C/j(n) in order to get a limit. Following 
|AP1| we assume that the total Hamiltonian, describing one elementary interaction, is of the 
form 

Htot = H®I + I(^ (^^ +^{C (^a!i + C* ®al) . 

That is, a typical dipole-type interaction Hamitonian with a renormalization in ^/n of the 
field operator and aj in order to strengthen the force of the interaction while the time of 
interaction decreases. 

With this Hamiltonian, it is easy to check that the coefficients of U{n) are of the form 

U^{n) = I + l(-iH-i-tQl+\c*c]+o{-) (10) 
n \ Z J n 

C/?(n) = -i^C + o(i) (11) 
Jn n 



Ul{n) = -i^C*+o(i) (12) 

(13) 



Ulin) = I + -(-iH-i-iiI+\cC*]+o{-). 

n \ Z J n 



3 Continuous Trajectories 

In this section, we implement the asymptotic expression of the coefficient Uj{n) in the descrip- 
tion of the quantum repeated measurements for the model of the two level atom. First we 
recall the convergence of discrete models to continuous models of Belavkin equations. Second, 
we show return to equilibrium results in this context. 

As in shown in |Pel| and |Pe2] , the continuous limit of the evolution equation is completely 
different, depending on wether the observable A is diagonal or not in the basis of rj. The point 
is that the limit equation is of diffusive type when A is non-diagonal and of Poisson type in 
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the diagonal case. Inside each case, the behaviors are very comparable and differ only by 
some coefficients. This is why, it is enough here to consider only two cases: 



A 

as representing the diagonal case, or 




1 



1 

1 



A=[ " ) = + a? , 



as representing the non-diagonal case. Here, we focus on the description of quantum trajec- 
tories in terms of pure states, while in |PeH IPe2] . the evolution for the density matrices is 
considered. 



3.1 The Poisson case 

We first start with the case A = a\, for which we have Pq = a^. It is easy to see that we 
can choose = 1, = for the description of the projectors Pi. Applying the hypothesis 
(fT0]l - (fT3]l . we obtain the probabilities 

Pk+i = Tr{pkPa] = \\UQ\'4)k)\\ = l--]:lik{n)+o 

n I 

Qk+i = Tr[pkPi] = \\U^\il^k)\\ = - 7;l^k{n) +o 

n 2 

where Hkin) = {ipk, C*C ipk)- By remarking that 1q = 1 — Ij^, we have the following difference 
equation for (ipk)'- 

\tPk+i) - m = \ (^-iH - ^C*C + i/ifc + o(l)^ 

+ (-^-i+o{i))mi',^'. (14) 

In the continuous limit, we shall see that this difference equation converges to an equation 
of the form 

d\^l^t) = [-iH - ^ {C*C + fit- I) + VJ^C^ \i^t-) dt+ 

+ - X- \i^t-) {dNt - fit- dt) (15) 

where pt = {i^t,C*C tpt) and (Nt) is a counting process such that t — )• iV^ — f^fXgds is a 
martingale. This is to say that (Nt) is a counting process with stochastic intensity equal 
to Jq* fig ds. A first problem is that equation (fT5]) is ill-defined. Indeed, the intensity of the 
counting process depends on the solution itself. We need to be more precise about what we 
mean by a "solution to equation (fT5|) ". 
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Definition 1 Let (fi, J^, P) be a probability space. A process-solution of the jump-equation 
([7^ is a process (V't) and a counting process Nt, with intensity fig ds where fit = {^Pt, C*C tpt), 
such that for all t we have 

\^t) = li'o) + (^-iH - ^ {C*C + fis-I) + VJITC^ \^Ps-)ds+ 

+ r ^^~'^^\ i^s-){dNs-fis-ds). (16) 

This notion of solution imposes the simultaneous existence of the process and the 
counting process Nt- In order to construct such a counting process, we use a Poisson point 
process. 

Let {Q,J-,P) be a probability space, on which is living a Poisson point process N on 
such that the expectation of the number of points N{u},B) lying inside a Borel set B is given 

by 

E[N{-,B)] = X{B) 

where A is the Lebesgue measure on M?. 

This way, N defines a random measure B i— )• N{(jj,B) on M?, whose volume element is 
denoted by N{uj, dxds). The following theorem shows how the random Poisson measure is 
used to construct the counting process. 

Theorem 5 ( |Pe2j ) Let (17, J^, 7"^, P) be a filtered probability space on which lives a Poisson 
point process N . The following equation 

li^t) = IV-o) + 1^ [-iH - \ {C*C - fis- /)) IV's-) ds+ 

+ f I - ^gp-^) |^^_^ io<,<^^_ N{dx, ds) . (17) 

Jo Jm. v^«- 

admits a unique solution (tpt) such that WiptW = 1 almost surely. Furthermore the process 
(l^i)) together with the counting process 

Nt= [ f lo<.<^.,_N{dx,ds) (18) 
Jo Jn 

constitute a process-solution for equation US]). 

Even if this theorem is just an application of the results of [Pe2j . let us explain roughly 
how it is proved (this description will allow also to describe the return to equilibrium property 
in the jump case). 

In equation (|17|) there are two parts: the ordinary differential part and the one driven by 
the Poisson process. Consider the collection of jumping times of the Poisson process. If there 
is no jump of the Poisson process A'^, we deal with an ordinary differential equation 

li^t) = IV-o) + 1^ [-iH - \ {C*C - fis- /)) \i's~) ds . 
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This equation admits a unique solution, from which we deduce the curve t ^ pt- The first 
time Ti when the Poisson process has a jump under this curve, the solution jumps and 
takes the value 

After this first jump, we have a new "initial" value for \^/Jt) and the process starts again in 
the same way: we solve the ordinary differential equation and the solution follows it, until it 
meets a jump of N which is bellow the curve, then it jumps. And so on. 

Remark: The corresponding evolution for the density matrices can be obtained by com- 
puting the stochastic differential equation for pt = \tpt){'>Pt\- By applying the stochastic 
calculus rules for random Poisson measure, we get the equation 

Pt = Po + {l{ps-) - Cps-C* + Tr[C/5,_C*]p,_) ds+ 

where L is the Lindblad operator defined by 

L{p) = -i[H,p]-]^{C^C,p] + CpC\ 

Thus, by defining J{p) = CpC*,we recover the equation ([1]) mentionned in Introduction. 

Now that equation (jl7p is well understood, we wish to pass to the continuous time limit 
on equation ()15p . The appropriate topology for the convergence theorem proved in |Pe2] is 
the Skorohod topology. Let us recall it. For all T > we denote by P([0,T]) the space of all 
cadlag matricial process on [0, T] endowed with the Skorohod topology, that is, the topology 
of the weak convergence of cadlag processes (the convergence in distribution). 

The approximation result is based on the description of a quantum trajectory as the 
solution of a stochastic equation wich is a discretization of (jl7p . In particular, from equation 
we can write 

[nt]-l 

iv'm) = i^o>+ Yl (i^fc+i>-i^fc») 

A:=0 

= i^o) + -i-iH- -c*c + -pk + o(i) J m + 

+ f-^-/+°(l))lV'.)l^^ (20) 

for all t > 0. An adaptation of the result of |Pe2j give us the following convergence. 

Theorem 6 ([Pe2j) Let T be fixed. Let {Q,J^,P) be a probability space in which lives a 
Poisson point process N. Let {\'4'[nt]))o<t<T be the discrete quantum trajectory defined by 
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the equation l\20\} . This discrete quantum trajectory converges in P([0,T]) to the process 
i\ipt))o<t<T which is the unique solution of the stochastic differential equation 



C 



where fit = {tpt , C*Ct[)t). 

3.2 The diffusive case 

We now consider the case where 



A 



)lo<x<Ms- N{dx,ds) 
(21) 




'1 1^ 

2 2 



1 1 
v2 2/ 



1 1 



We have 



Pn 



1 1 

2 2 



1 1 

2 2 



and fi = u = Hence, after computation we obtain: 

Pk+i = Tr[pkPo] 
Qk+i = Tr[pkPi] 



i=(C/0-[/0)|V.) 



1 




+ o| 




2 + 






y/n 




Vn 


1 




+ o| 








2 ~ 


\fn 




V n 



where ffc(n) = Re(V'fc, CV'A;)- 

We introduce the random variables [X]^ defined by 



for all A; > 0. In terms of the evolution equation takes the form 



1^, 



IV'fc) = ^ i-iH - ^{C*C - 2ukC + ull) + o(l) ) 1^^^) + 



(22) 
(23) 



k+l) 



The continuous diffusive equation which is the natural candidate to be the limit of equation 
(IMD is 



C-uu + o{l)]m—Xk+i. (24) 
n 



d\tlJt 



{-iH - 1 [C*C - 2utC + lyfl)^ \^Pt)dt + (C - i^tmt)dWt, 



(25) 



where ut = Yie{tpt,Cipt) and {Wt)t is a one-dimensional Brownian motion. 

In [Pel], it is shown that the convergence result is highly based on the existence and 
uniqueness of the solution for such equation (let us stress that the coefficients are not Lips- 
chitz). In particular, by a truncation method the following Theorem is proved in |Pel| . 
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Theorem 7 ( |Pelj ) Let {n,J^,J^t,P) be a probability space on which is defined a standard 
Brownian motion {Wt)t- The following stochastic differential equation 

d\i^t) = {^-iH - ^ {C*C - 2vtC + vp)^ \^t)dt + {C - utl)\i^t) dWt (26) 

admits a unique solution. Furthermore, almost surely, for all t we have \\il^t\\ = 1- 

We can now consider the approximation procedure. In a similar way as the Poisson case, 
we can consider the difference equation 

[nt\-l 

fc=0 

[nt\-l 



E l^{-iH-\[C*C-2vkC + vll) + o{l)\\i,^) 

[nt]-l 

+ (c\i^k)-'^km+oii))^Xk+i. (27) 



k=0 

We have the following result. 

Theorem 8 Let T be fixed. Let {Q.,T ,J-t,P) be a probability space on which is defined a 
standard Brownian motion {Wt)t- Let {\'4'[nt]))o<t<T be the discrete quantum trajectory defined 
by the equation (27). This discrete quantum trajectory converges in P([0, T]) for all T to the 
process (|V't))o<t<T which is the unique solution on O of the following stochastic differential 
equation: 



(28) 



d\^t) = {C- vtL)\^t) dWt + {-iH - i {C*C - 2vtC + u^j)^ \^t) dt 
where vt = ^^ifpt^CTpt) . 

For sake of completeness we give some details on how to prove such a convergence. 
Proof: Define the processes 

Mt) = IV'm), Vn{t) = Wn{t) = ^ V Xfc+i. 

The process {ipn{t)) satisfies 

Mt) = (^-^C*CMs-) + Re{i'n{s-),CMs-)))CMs-)^ dVnis) 

+ l\ci>n{s-) - ReiMs-), CMs-Wn{s-) dWn{s) + Sn{t), (29) 







where the terms en{t) corresponds of the o terms in the equation in asymptotic form. 
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In order to show that V'n(i) converges in the Skorohod space to the solution of 

IV-t) = 1^ (^-^C*C\iPs) + lytCli^s)^ ds + l\c\4^t) - i^M) dWs 

we make use of the celebrated Kurtz-Protter theorem. Let us recall it. 

Recall that [X, X] is defined for a semi-martingale by the formula [X, X]t = X^ — 
Jq Xg- dXg- For a finite variation process V we put TtiV) to be the total variation of V 
on [0,t]. 

Theorem 9 (Kurtz-Protter, [K-P| ) Suppose that Wn is a martingale and Vn is a finite 
variation process. Assume that for each t > O." 

SUpE[[Wn,Wn]t] < (X) 
n 

supE[rt(K)] < oo 

n 

and that {Wn^Vm^n) converges in distribution to {W,V,Q) where W is a standard hrownian 
motion and V{t) = t for all t. Let Xn{t) be a process satisfying 

Xn{t) = po + Snit) + f L{Xn{s-)dVn{s)+ f Q{Xn{s-)) dWn{s) 

Jo Jo 

Suppose that X satisfies: 

Xt = Xo+ I L{Xs)ds+ ['e{Xs)dWs 
Jo Jo 

and that the solution of this stochastic differential equation is unique. Then Xn converges in 
distribution to X . 

In our case, the different hypothesis above are satisfied. Indeed, define a filtration for the 
process {Wn{-))'. 

j;" = a(X„i< [nt]). 

The following is proved in [Pel]. 

Proposition 3 We have that J-""") is a martingale. The process (Wn{.)) converges to 

a standard Brownian motion W, when n goes to infinity and supn'E[[Wn,Wn]t] < oo. 

Furthermore, we have the convergence in distribution for the process (Wn, Vn, £n) to (W, V, 0) 
when n goes to infinity. 

This proves the announced convergence. □ 
Remark: Using Ito rules on \'iljt){ipt\, we get the equation for density matrices 



dpt = L{pt) dt + [Cpt + ptC* - Tr[(C + C'')pt]pt) dWt, 

which corresponds to the equation ([2|) mentionned in the Introduction (the Lindblad operator 
L has the same form as the Poisson case). 
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4 Return to Equilibrium 

Now that the hmit equations are estabhshed, we are interested into the long time behaviour 



of the solutions. We specify our investigations to the special case where C 
1 







1, , 



^0 



and H = Hr i q q 

Writing the processes (^f) in terms of their coordinates; that is (V't := {xt,yt)), the 
Belavkin equations take the form 



Xt = Xo + 



j i^-ixs + Re{xsys)ys - ^Re{xsys)x^ ^ J ~ ^^^^^y^^^''^ 
yt = yo + J (^-^Vs -^Re{xsysfys^ ^^^J (^~^^(^>^ys)ys^ dWs 

(30) 



in the diffusive case, and 
rt 



Xt 



xo+ I [ -ixs + \xs\ys\^ \ ds+ [ [ i - Xs- + I] N{dx, ds) 

Jo \ ^ J Jo Jo<x<\y,^\2 \ J 

(31) 

yt = yo+ I (-T^ys + hysl'^ys] ds+ I [ i - ys-] N{dx,ds) 

Jo \ ^ ^ J Jo i0<a:<|j/,_|2 \ ) 

in the Poisson case. 

In the Poisson case, note that the intensity is lit = \yt- P, so that one can restrict ourselves 
to the case where the jumps of the Poisson process are in between the lines y = 1 and y = 
(we have namely |yt-P < 1, for all t). The function t — )• card{N{., [0, 1] x [0, t])) = Mt then 
defines a standard Poisson process with intensity 1. The Poisson random measure and the 
previous process generate on [0, T] (for a fixed T) a sequence {(Tj, ^j), i G {1, . . . , A/i)}} where 
each Tj represents the jump time of A^. Moreover the random variables are uniform random 
variables on [0, 1]. Consequently we can write our quantum trajectory as follows 



a^o + ^ y-ixs + ]^Xs\ys?j ds + Y^^y- XT,- + 1 j lo<c,<|yT,-| 
yt = y^ + + \\ys?y^ ds + Y,(^- yx^^j lo<^^<|j,^^_| 



(32) 



Now, we shall investigate the large time behaviour of a solution of equation ()30p or ()3ip . 
To this end we need to notice the following lemma. 

Lemma 1 Let ^ "^^ ^ 6e the either solution of equation 130)) or \31\) starting with an initial 
condition of the form f j • Then, almost surely, we have yt = 0, for all t. 
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Proof: Starting from yo = 0, in each case, it is easy to verify that yt = is a particular 
solution for the corresponding stochastic differential equation describing the evolution of {yt)- 
As a consequence by uniqueness of solution, almost surely, yt = for alH. □ 

Remark: In both cases, if yt = for all t, it is easy to see that the evolution of {xt) is 
given by the solution of dxt = —ixt dt. 

Remark: In terms of states, this lemma expresses that if ■0o = ( J , we have almost 



surely 







for all t. In other words, the state is an invariant (or stationnary state) for the 

stochastic dynamic of continuous measurement (let us stress that without measurement, i.e 
in the deterministic regime, it is easy to sec that this state is already the invariant state). 

Now we can make precise the result which states the return to equilibrium property. In 
particular we focus on the large time behaviour of the part yt and we show that this process 
converges to zero when t goes to infinity. 

Proposition 4 Let \il)t) = ( ^* ) 6e either the solution of the jump- equation or the solution 



. yt _ 

of the diffusive equation, then we have 



Therefore, we have 



\ytV^Q. (33) 



(34) 

t— >-oo 



and the process of pure states (|V't)(V't|); where i^t = \ ^* ], for all t, satisfies 

\ yt J 

|Vt)(VtL-^|f^)(f^|- (35) 



Proof: Let us first treat the case of the jump-equation. We need to share into two cases, 
if there is jumps or if there are no jumps. 

In the case where there is at least one jump. At the first jumping time Ti we have 

XTi \ _ f XTi- \ _^ f -XTi- + \ _ f 1 



Following the description of the solution of the jump equation, the solution after Ti is given 

by the ordinary differential part with the new initial condition ^ q ^ ■ This initial condition 

satisfies yo = 0, then by Lemma 1, we get yt = for all t > Ti. 

If there are no jumps, this corresponds to the event A = {to ^ ^/N{lo, {{s,x) G < 
X < |j/sP} = 0}. In this situation, the evolution of (yt) is only given by the ordinary differential 
equation 



1 r 

2/t = yo + 2 i-ys + bspys) ds . 
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We want to show that \yt\^ — )• 0, when t goes to infinity. Derivating, we get 

By Lemma 1, if yo ^ 0, we have \yt\'^ > 0, for all t. In this case, dividing by y| we solve the 
equation and we get 



= ly^P X exp [ -{t- s)+ / \yufdu 



for all t > s. In particular the function t — )• is decreasing, then we get 

|yt|2 < |y,|2exp(-2(t-s) + 2(t-s)|y,|2) , (36) 

for all t < s. Since we have l^^p + ly^P = li for all s, we have \ys\'^ < 1, for all s. 

With the estimate (p6]) . in order to conclude, we need to show that there exists s such 
that ly^p < 1. Assume the contrary, we should have \ys\ = 1, for all s. But such an event is 
of probability zero. Let us prove this fact. In this case the event A becomes 

A = {uj €n; N{lo, {{s, x) G < x < 1, s > 0}) = 0} . 

We claim that A is of probability 0. Indeed, we have 

P[{uj e n ; N{uj, {(s, x) G M^^ < X < 1, s > 0}) = 0] = 

= lim P[{uj G Q ; N{uj,{{s,x) € < x < 1, < s < n}) = 0}] (37) 

n 

and 

P[{uj G n ; N{lo, {(s, x) G I^VO < x < 1, < s < n}) = 0}] = exp(— n) . 

Hence the announced result about A. 

As a consequence, there exist s such that ly^P < 1- For this s, by taking the limit t goes 
to infinity in expression ([36|) . we get \yt\'^ — )• 0. 

With the above discussion, for the jump equation, it is easy to conclude that 

2 a.s 

— ^ . 

I— >CXD 

Let us now treat the diffusive case. In order to prove the result we shall show first that 
converges almost surely to a random variable Uoa when t goes to infinity. Second we 
show Uoo = almost surely. Using Ito rules, we get 

d\yt\'^ = ytdyt + yt dyt + dyt dyt 

= -\yt\'^dt-2Re{xtyt)\yt?dWt 

As a consequence we have almost surely: 

yl = yl+ [ -\yuf du+ f -2Re{x-uyu)\yuf dWu, (38) 

J s J s 
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for all t > s. Let (J't) be the filtration generated by the Brownian motion, that is J-f = 
a{Wu,u < t}. Since E —2Iie{xuyu)\yu\'^ dWu\J^s = 0, the above equation shows that 

E[\yt\^\Ts]<n\ys\^]- 

This way the process is a super martingale which is bounded (for all t, we have < 

\yt\ < !)• Therefore, this process converges almost surely to a non- negative random variable 
«oo when t goes to infinity. In order to show that this random variable is equal to zero almost 
surely, we just have to show that E[uoo] = 0. To this end, from Eq. [38] for s = 0, we get 

n\yt\^]=yo+ f-mysf]ds. 

Jo 

Solving the equation, we get 

E[|yi|'] = lyol'e-*. 

As a consequence, we get 

E[y?]^^ 0. 

Now, using the Lebesque dominated convergence Theorem, we deduce that E[moo] = and 
then = almost surely. The proposition is then proved. □ 

Remark: In the proof, we have supposed that the initial condition is deterministic. This 
result can be easily genralized by assuming that the initial condition is random and the same 
result holds. 



Remark: In Probability Theory, usually we consider invariant measure for stochastic 
process. Here the invariant measure is the Dirac measure on the state In |BaPj . es- 

sentially for the diffusive equations, by assuming special conditions on the coefficients defining 
the stochastic differential equations, results for invariant measure for stochastic Schrodinger 
equations have been investigated. They established results where the invariant measure owns 
particular properties (in particular concerning the support of the measure which is absolutely 
continuous with respect to the Lebesgue measure). 
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